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We argue, using the AdS/CFT correspondence, that the transient dynamics of the shear stress 
tensor in a strongly coupled Af = 4 SYM plasma is not described by relaxation-type, fluid dynamical 
equations: at long times the equations of motion should contain a second-order comoving derivative 
of the shear stress tensor. This occurs because in this strongly-coupled system the lowest "non- 
hydrodynamical" quasinormal modes associated with shear stress possess a nonzero real part at zero 
wavenumber. We use Weyl invariance to obtain the most general equations of motion containing 2 
comoving derivatives of the shear stress tensor in the transient regime that are compatible with the 
symmetries. We show that the asymptotic solution of this theory valid at times much larger than 
the timescale associated with the "non-hydrodynamical" modes reproduces the well-known results 
previously obtained directly from the AdS/CFT correspondence. If the QGP formed in heavy ion 
collisions can be at least qualitatively understood in terms of strongly-coupled J\f — 4 SYM theory, 
the second time derivative present in the equations of motion of the fluid may lead to an unexpected 
dependence on the initial conditions for the shear stress tensor needed in numerical hydrodynamic 
simulations. 
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I. INTRODUCTION 

Relativistic fluid dynamical models have played a key role in our current understanding of the nearly perfect fluid 
behavior displayed by the Quark-Gluon Plasma (QGP) formed in heavy ion collisions In these models, exact 
energy-momentum conservation, V^T^ — 0, is supplied by another phenomenological dynamical equation for the 
macroscopic shear stress tensor, 7r Miy , which is defined as 

n J *' = r^-EtiY + PA' B ' 1 (1) 

where T^ v = (T MI/ ), e is the local energy density, P is the local pressure, u M is the local fluid 4-velocity, and 
A^ v = g^ u — u^u v is a spatial projector (our metric signature is +, — , — , — ), i.e., u^A^ = (we shall not consider 
here the contribution from the bulk viscous pressure or the effects from nonzero baryonic density). Since we are 
going to consider only conformal fluids in this paper, the trace of tt^ is equal to zero. In fact, in terms of the 
doubly symmetric and traceless projection operator, A^ Q/5 = (A^ Q A^ + A^A" a ) /2 - A^A afl /3, one can see that 
A^tt"/ 3 = tt^ and u^ v = 0. 

In relativistic fluids causality is intimately connected to stability @, Q and Israel and Stewart [3| were among the 
first to understand that the characteristic times within which fluid dynamical dissipative currents, such as ■n liV , relax 
towards their asymptotic Navier-Stokes values cannot be arbitrarily small. The so-called Israel- Stewart (IS) equations 
for the shear stress tensor are relaxation-type equations (in terms of the co-moving derivative d/dr = w^V^) of the 
following general form 

n jr<""> + = 2r]<T f " / + . . . , (2) 

where t\ is the relaxation time coefficient, rj is the shear viscosity, a^ v = V^u^ is the shear tensor, = A^V 1 ', 
A^ v ^ — A^dA a ^/dr, and the dots denote nonlinear terms involving tt^ and gradients of T and [3, The 
major step taken by Israel and Stewart in Q was to realize that causality demands that the dissipative currents 
obey dynamical equations of motion (which introduce the linear transport coefficient t\) that describe their transient 
dynamics towards their respective asymptotic relativistic Navier-Stokes solution. A few years earlier, Kadanoff and 
Martin [6] argued that a similar relaxation transport coefficient should appear in the description of spin diffusion in 
a way that is consistent with well-known sum rules. 

Using the 14-moments method, it is possible to show 0, 0] that in relativistic gases t\ is of the order of the 
microscopic collision time. Therefore, one expects that in physical situations where the time variation of the fluid flow 
is comparable to this microscopic scale, the fluid is in the transient regime where the relaxation dynamics described 
by ti becomes important. For dilute gases, such as air, under normal circumstances the collision time is orders of 
magnitude smaller than the typical time variation of the flow. However, it is possible to create physical systems under 
which the flow of a fluid varies in a timescale of the order of the mean free time (such as in microflows 7] ) . Given the 
rapid expansion experienced by the QGP formed in heavy ion collisions, it is reasonable to investigate the dependence 
of hydrodynamic predictions on the actual value of T\ (see, for example, |8(). 

If the dynamical properties of the strongly-coupled QGP can be (at least qualitatively) understood using M = 4 
Supersymmetric Yang Mills (SYM) theory, we shall see in this paper that this would imply that the relaxation 
equations for ir^ commonly used in numerical simulations must be replaced by new equations for transient dynamics 
involving second-order comoving derivatives of . 

II. NON-HYDRODYNAMIC POLES AND TRANSIENT FLUID DYNAMICS 

It is well-known that retarded correlators can have singularities such as simple poles and also branch cuts. Of 
particular relevance for fluid dynamics are the so-called "hydrodynamic" poles, wo(k), which appear in retarded 
correlators of conserved currents Q. The k dependence of these modes can be used to obtain the corresponding 
diffusion transport coefficient D associated with a given conserved quantity through the relation lim^o s ) ~ 
— iDk 2 + . . . [9]. The hydrodynamic modes are characterized by the momentum dependence ~ —ik 2 and, consequently, 
vanish in the limit k — > 0. Because of their appearance in Navier-Stokes theory, the existence of hydrodynamic modes 
is quite often taken as evidence for fluid dynamical behavior. Furthermore, modes that do not share this behavior, 
i.e., modes in which limk-s-o w n (k) ^ 0, are known as "non-hydrodynamic modes". 

Surprisingly enough, a stable theory of relativistic fluid dynamics cannot be formulated using only hydrodynamic 
modes. In order to obtain a causal and stable theory of fluid dynamics, the shear stress tensor has to be promoted to 
an independent dynamical variable, as in Israel and Stewart theory, e.g. Eq. Since the shear stress tensor is not a 
conserved quantity, when it becomes an independent dynamical variable non-hydrodynamic modes must appear in the 
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theory. In the case of Israel-Stewart theory, the non-hydrodynamic modes describe the relaxation of the dissipative 
currents towards their respective Navier-Stokes asymptotic solutions and they can be directly related to the relaxation 
times (iq|. 

Clearly, in order to fully describe the complicated many-body dynamics of an interacting system down to arbitrarily 
small time scales, all the infinite number of non-hydrodynamic modes must be taken into account. For sufficiently 
long times, however, only the slowest modes should contribute significantly and one should be able to systematically 
neglect the effect of faster modes to the system's dynamics (at infinite times, no non-hydrodynamic mode should be 
required at all). This type of truncation should be possible at sufficiently long times and does not depend on whether 
or not the non-hydrodynamic poles are parametrically separated (as long as the distribution of poles is discrete) . 

As a matter of fact, it was shown in [10] that the long-distance, long-time linearized dynamics of the shear stress 
tensor in any system that can be described via the Boltzmann equation must follow the general ansatz from Israel 
and Stewart, Eq. ([2]). In other words, at long times ir^" must obey a first-order differential equation in the comoving 
derivative that describes how it relaxes towards its steady state Navier-Stokes value. This means that in relativistic 
dilute gases the relaxation time can be extracted directly from the microscopic theory and does not necessarily have 
to be considered as a regulator for the gradient expansion, as advocated in (for a recent discussion see [12(]). We 
remark that this relaxation behavior at long times should not be taken for granted: it is a direct consequence of the 
fact that the retarded Green's function associated with shear stress tensor, obtained via the linearized Boltzmann 
equation, is a meromorphic function where all the poles lie on the (negative) imaginary axis and, in this case, for times 
tu>i(Q) 3> 1 one obtains that T\ = —i/u>i(0) [Iol |. with u>i being the non-hydrodynamic mode with smallest frequency. 

The near-equilibrium dynamics of weakly-coupled QCD at very large T is expected to be described by a Boltzmann 
equation involving quark and gluon scattering |13| . Therefore, according to the general discussion in [lfj , one should 
expect that deep in the deconfined phase long distance, long time disturbances of the shear stress tensor follow IS- 
type equations although the correct value of t\ in this case, as determined by the first non-hydrodynamic pole of the 
retarded Green's function, has not yet been computed. 

Sufficiently below the phase transition, say T ~ 130 MeV, lattice data 14] for QCD thermodynamics seems to be 
consistent with the predictions from non-interacting hadron resonance models. Thus, it is natural to assume that at 
low T QCD behaves as a weakly coupled gas of hadrons and resonances, which then implies that a description via 
the Boltzmann equation (including different particle species) may be appropriate. Using the formalism derived in 
[lfj, it is possible to show that the total shear stress tensor in this case, which is a sum over all the hadronic species, 
follows an IS-like equation of motion with a single relaxation time coefficient given by the slowest non-hydrodynamic 
mode. Therefore, we expect that in QCD the non-hydrodynamic pole of the retarded Green's function closest to the 
origin should be a purely imaginary number at very high or low T and in the transient regime one will find relaxation 
equations for n^". 

However, it is clear that the surprising "perfect fluid" character of the QGP appears not at low or high temperatures 
(where r//s ~ 1) but rather near the phase transition, say T ~ 160 — 200 MeV, where the number of degrees of freedom 
increases very rapidly and r//s, the figure of merit for perfect fluid behavior, is expected to be rather small [15j . Given 
that the correct physical mechanism that leads to the perfect fluid behavior of the QGP is not yet known and that 
the relevant t'Hooft coupling near the phase transition, Xqcd = 9qcd ~ 10 for N c = 3 and olqcd ~ 0.3, we find 
it useful to investigate what the AdS/CFT correspondence has to say about the fluid dynamical equations for 
the shear stress tensor in strongly-coupled gauge theories in the transient regime. In other words, how do strongly 
coupled plasmas (which possess gravity duals) relax towards their asymptotic, universal Navier-Stokes solution? This 
question will be answered in the following sections. 



III. WEYL INVARIANCE AND THE EQUATIONS OF MOTION FOR THE M = 4 SYM FLUID IN THE 

TRANSIENT REGIME 

Recently, a new way to derive the equations of motion of relativistic fluid dynamics based on Weyl invariance was 
put forward by Baier et al in Ref. [TTI j . The main idea is to use the fact that the dynamics of conformal plasmas (with 
equations of motion involving less than 4 derivatives) should be invariant under Weyl transformations in which the 
metric changes as g^ LU — > g lliy (x) e~ 2n ( x ) ; and Q(x) is an arbitrary scalar function. Since the energy momentum tensor 
scales classically, it is easy to prove that under a Weyl transformation it scales homogeneously with conformal weight 
equals 6, i.e, T^ v — > e 6 °T A " y . The basic hydrodynamic variables change under Weyl transformations as follows: from 
u^u^ — 1 one obtains that — > e n and, since the ideal energy-momentum tensor has conformal weight equals 6, 
the temperature scales as T — > e a T and, thus, the dissipative part of the energy-momentum tensor also transforms 
homogeneously, i.e, — > e 6n TT^. 

It is important to notice, however, that while T and scale homogeneously, their conventional spacetime covariant 
derivative does not. Things get significantly easier with the aid of the Weyl covariant derivative defined in [l7j . which 
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acts on the basic hydrodynamic variables as follows 

u a 9 T 

D a T = V a T — T u a H g — , 

A V Q 

D a u u = V a u u - u a u' / ^- , (3) 

where 9 = V ' v u v . In this case, these derivatives are homogeneous under Weyl transformations, i.e., D a T — > e n D a T 
and D a u v — > e n D a u v . Defining the comoving Weyl invariant derivative as D = u^D^, one can see that Du^ = 0, 
D^u^ 1 — 0, and DT = T + T9/3. Moreover, since tt^ is traceless and of conformal weight equals 6, it is possible to 
show that D a n a P = V Q 7r Q ^, A£V a ir alJ = D± a TT al3 + u^Tr a p(j a P and we can now write the general conservation laws 
for a conformal fluid (where e = ip) as 

3 ' 

D^P = D ±a ir a » + u^n aP a a ^ (4) 
where DP = P + 4P6»/3 and D^P = V^P - 4?^. 



A. The Gradient Expansion Approach in Conformal Fluid Dynamics 

In order to solve the conservation laws (|4]) we must provide the equation satisfied by the shear stress tensor. One 
possibility to derive such equation is via the gradient expansion in which is assumed to be solely expressed in 
terms of P (or temperature), and their gradients. In this framework, it is possible to express in terms of a 
controlled expansion in powers of derivatives or order of derivatives of P and , 

tt^ = m irf + % nf + • • • , (5) 

where the quantities TL^" and IT^ correspond to terms of first and second order in gradients of P and it M , respectively, 
and the dots denote possible terms with higher order derivatives. 

This derivative expansion is controlled by a small parameter called the Knudsen number, Kn = ^micro/^macro, 
which is basically the ratio between a microscopic length scale (e.g., the inverse temperature for conformal fluids or 
the mean free-path for gases) <~ ^ m icro and the overall macroscopic length scale of the fluid ~ I macro (the inverse of 
the gradient of velocity or temperature). The term IK J " y is assumed proportional to the gradient of a macroscopic 
variable and should be of order ~ -^macro- Every additional derivative brings in another inverse power of L maC ro and, 
thus, TL^ V ~ L~" cro . The microscopic scale £ m i cro is contained in the coefficients r\i. Up to some overall power of 
^micro (which restores the correct scaling dimension), rj n ~ ^micro- Therefore, the terms IK J " y and multiplied by 
their corresponding coefficients in Eq. ([5]) are of order Kn and Kn 2 , respectively. Subsequent terms would be of a 
higher order in Kn and, therefore, when the system exhibits a clear separation between Anicro and L maC ro, he., when 
Kn <C 1, it is possible to truncate this expansion. Ideal fluid dynamics corresponds to the zeroth order truncation of 
this series, i.e., when no terms are included at all. 

The first term in the gradient expansion can be obtained by constructing all possible tensors that can be made 
using first order derivatives of P and it M . These can be easily obtained and are: 

D^P and D^u u . (6) 

Next, one has to build, using these gradients, tensors that have the same properties satisfied by and, at the same 
time, transform homogeneously under Weyl transformations. Here, the only possibility is the shear tensor, which 
can be written in terms of the Weyl derivative of u M as 

a^=A^VX= \ ■ (7) 

Therefore, the most general equation allowed by symmetry that can be satisfied by 7r MI/ , up to first order in Kn, is 

which corresponds to the relativistic Navier-Stokes theory, with r)\/2 being identified as the shear viscosity coefficient 
■q. It is easy to see that the shear tensor scales and, therefore, i] ~ T 3 . 
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In the framework of the gradient expansion, relativistic Navier-Stokes theory can be extended by including terms 
of second order in gradients of P and u M . In order to do so, one has to obtain all the possible terms involving Weyl 
derivatives of P and that contribute to II2 " ■ All the independent terms of second order in gradients of pressure P 
and w M that are symmetric, transverse, traceless, and that transform homogeneously under Weyl transformations are, 

Of = Da^ v) = + (7^9/3, 
Of = P> Il/) + 2u a R a ^ l3 u , 

of = o-j[V> A , of = a { fn v », of = n<$p?> x , (8) 

where ft^ = (V^w" — V±u fM )/2 is the vorticity operator, R^" is the Ricci tensor, and R^ vot P is the Riemann tensor. 



All the terms above have conformal weight 4 and were first found and listed in Ref. Note that terms such as 



A^DfD^P, A^(DfP)(D^_P), and A^(DP)a a/3 contribute only to 0(Kn 3 ), as can be seen by substituting the 
leading order relation tt^ ~ a^ v together with the general conservation laws (QJ. 

Therefore, the most general equation allowed by symmetry that can be satisfied by 7r A " y , up to second order in Kn, 

is 



2na» v - 2r l b * °i V + ° ( Kn3 ) • ( 9 ) 



The 6 coefficients r\ and hi can be calculated via Kubo formulas for the correlators of the energy-momentum tensor 
derived using metric perturbations. In strongly-coupled M = 4 SYM theory, all the coefficients above (including those 
associated with nonlinear terms) were determined using the AdS/CFT correspondence [U El Hjjj . For instance, for 
strongly-coupled SYM one finds t](2ttT) / Pq — 2 (Pq ~ T 4 is the pressure at equilibrium) and &i(27rT) = 2 — In 2 [llj . 



B. Going Beyond the Gradient Expansion via the Inclusion of Transient Effects in Relativistic Fluid 

Dynamics 

Relativistic Navier-Stokes theory and its extensions via the gradient expansion are hindered by acausal behavior 
which complicates their usage in relativistic problems [2j . In Ref. , a stable and causal fluid-dynamical theory was 
obtained from the gradient expansion by substituting in all second-order terms the "inverted" first-order solution, 
a^ u — n**" / (2rj) . Then, the following equation of motion for tt^ appears, 

biDic^ + 7T^ y = 2r/a^ - 2r)b 2 Of - 27763 Of 

- 27764 Of - 2 V b 5 Of , (10) 

where Of corresponds to Of with the substitution — > i:^" /(2r)) and we used that DT ~ 0(Kn 2 ). Note, 
however, that in order to render the gradient expansion stable, the shear stress tensor had to be promoted to an 
independent dynamical variable. On the other hand, Eq. Q was proved to be the most general equation allowed by 
symmetry only when ■n liV was not an independent dynamical variable. Therefore, for causal theories of fluid dynamics 
the analysis first proposed in Ref. [ll| and reproduced in the previous section has to be revisited. 

In this section, we use Weyl invariance to obtain the full set of nonlinear differential equations that describe a 
conformal fluid when the shear stress tensor is a dynamical variable. Now, the idea is to extend Navier-Stokes 
theory by including all possible terms that can be constructed from gradients of P, u^, and that are symmetric, 
transverse, traceless, and that transform homogeneously under Weyl transformations. Then, in addition to the terms 
constructed in the previous section, we can also build new terms, e.g., 

Dir^ v \ L>V^\--- (11) 

ir£o»> a , n^", ir£n"> . (12) 

Note, however, that the shear stress tensor can no longer be expressed in terms of a series in powers of Knudsen 
number. By including terms of the form, e.g. Dtp'"') and D 2 ^**") , the shear stress tensor satisfies a partial differential 
equation and its relation with the Knudsen number is dynamical, as happens with Israel- Stewart theory, and not 
algebraic, as occurred in the gradient expansion. We organize the most general equation of motion for tt^ in the 
following form, 

• • ■ + X 2D 2 TT^ ) + Xl DTT^) + 

5 

= 2 W ^ + eiTT^ a + e 2 7r^> a + eaTr^fi") a - £ 2 V a Of + • • • , (13) 
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where the dots denote additional possible terms. Note that, 

3 

D 2 TT^ = ~ 2u p TT P ^U^ + —TT^e 2 + 39 TT^ + -TT^ 9 . (14) 

9 3 

The truncation of Eq. (|T3"|) is not trivial, as was discussed in Ref. 10]. The terms on the right hand side serve 
as source terms for the shear stress tensor while the terms on the left hand side describe the relaxation/oscilation of 
the shear stress tensor when perturbed by gradients. The right hand side of Eq. (p~3|) . i.e., the source terms, can be 
organized as a series in Knudsen number and the so-called inverse "Reynolds number" Re -1 = In^n^] 1 / 2 /Pq [2C| . 
Since 7r MI/ is an independent dynamical variable, the inverse Reynolds number can be considered as an independent 
small parameter that gives additional information on how equilibrium is approached. Therefore, it is possible to 
systematically organize the source terms of the equation of motion as an expansion in powers of both Kn and Re -1 . 
In this case, the terms eiiTa'a") 01 and e^TT^Q^ a , the term e^'Ka^ 01 , and the terms r]Ci0^ u , are all the possible 
terms of order C(Re _1 Kn), 0(Re -2 ), and 0(Kn 2 ), respectively. If we wish to describe the source terms only up to 
order 0(Re -2 , Re _1 Kn, Kn 2 ), they are enough. 

The truncation of the left hand side is more complicated since it cannot be organized as an algebraic series in powers 
of small quantities, such as Knudsen number or inverse Reynolds number. However, the order of the differential 
equation in the comoving derivative on the left hand side is equal to the number of non-hydrodynamic modes included 
in the dynamical description of the system. For example, if we include only the first comoving derivative of ~k^ v , e.g. 
Dn^v) j we have only one non-hydrodynamic mode, while if we also include the second order comoving derivative we 
would have two non-hydrodynamic modes. 

The main purpose of the gradient expansion is to correct Navier-Stokes theory in cases where the Knudsen number 
is not very small, i.e., the microscopic scale is no longer very separated from the macroscopic scales of interest. By 
including second order gradients of P and u M , the Navier-Stokes theory is extended to describe the dynamics at larger 
wavenumbers or smaller wavelengths. However, if the separation between the microscopic and macroscopic scales is 
no longer optimal, it is not enough to extend the applicability of the theory to describe higher wavenumbers, but one 
should also extend it to describe higher frequencies. This is the role played by the left hand side of the Eq. (|13l) . 
When more comoving derivatives of ir^ are included, more non-hydrodynamic modes are introduced in the theory, 
and a description of higher frequencies is obtained. Therefore, the structure of the left hand of Eq. (TT3"|) has to be 
determined by carefully matching the modes introduced in the macroscopic theory with the modes of the underlying 
microscopic theory, taking into account what are the relevant frequencies in the macroscopic domain. Also, one has to 
make sure that such matching can be done, i.e., the modes included in the macroscopic theory exist in the microscopic 
one. 

For dilute gases described by the Boltzmann equation, all the non-hydrodynamic modes lie on the imaginary axis 
in the complex w-plane (in the limit of vanishing wavenumber) [To| . In the long-time limit it is only necessary to 
include the non-hydrodynamic mode with the smallest frequency (at zero wavenumber) and the equation of motion 
for i: pv with source terms up to order 0(Re~ 2 , Re _1 Kn, Kn 2 ) becomes 

5 

Ti Dir^ + tt" u = 2r]<j^ + e^a v) a + e 2 ^n u) a + e^tt^ a - ^ 2 VCl Of. (15) 

i=l 

One should remark that, in general, the coefficients Cj's in the equation above are different than the 6j's in Eq. ©. 
One can see that the transient theory defined in Eq. (|15[) reduces to the well known result © in the limit of vanishing 
relaxation time. In this limit, we can obtain an asymptotic solution for 7r Miy by substituting the first order solution 
tt^ ~ 2r)o plJ into all terms in Eq. (fig]), which then implies that, asymptotically, t v Dtt^ ~ 2r)T„ Da^ v) + 0(Kn 3 ). 
In fact, one can relate the new coefficients with those in Eq. ^ as follows: b\ = t\ + c\, b 2 = c 2 , 63 = C3 — e\ — 
64 = C4 — e 3 , and 65 = C5. Therefore, Eq. (|T5|) leads to the appropriate asymptotic limit up to C(Kn 2 ). Also, 
one can show that the general theory (in flat spacetime) obtained from the Boltzmann equation using the moments 
method, as recently derived in [20], has the exact same form as (fT5j) in the conformal limit (massless limit and 
cross section a ~ 1/T 2 ). Note also that Eq. (|10l) can be seen as a particular case of our transient theory in which 
e\ = c\ = C3 = C4 = 0. 

Equation (| 1 3[) can be extended by including one more comoving derivative of ir IJ-1 ' , 

X2 D 2 it^ + xi Dtt { ^ + = 2 W ^ + a + e 2 w^ a + e 3 7r<f fi^ a 

5 5 

- £ 2 VCl + fi rt&P " -tD^D^) \ (16) 

i=l i=l 
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where D^D x ir v) x is the only source term found of conformal weight 8, 

D^D X 7r u > x = (V Q - 6ii a ) V A 7r A/3 . (17) 

Above, we included source terms of order 0(Re _1 Kn), 0(Re- 2 ), C(Kn 2 ), and C(Re _1 Kn 2 ). In principle, we could 
have also included source terms of order 0(Re -3 ) and C(Kn 3 ), but this is time consuming and outside the purposes 
of this paper. Again, it is easy to see that the theory displayed above reproduces the C(Kn 2 ) gradient expansion 
obtained in [ll[, since all the new terms included are of order C(Kn 3 ) when the asymptotic solution ~ 2^cr A " y is 
substituted. 

It is interesting to observe that, since the transient theories in (fl"5|) and (fT6|) reduce to (J9j) , several results previously 
derived using ([9]) are automatically valid also for the transient theories derived here. For instance, the expansion 
around k — > for the sound mode present in the theories defined via Eqs. (IT51) . (fTI)]) , and © is 



k 



«— (t) - ^ " '*bs ± £s - + ° [kt) ' (18) 

As mentioned above, the transient theory obtained from the Boltzmann equation taking into account only the 
slowest non-hydrodynamical mode assumes the form of Eq. (fl5|) . We shall see in the following that for the strongly 
coupled SYM fluid the transient theory assumes the form displayed in Eq. (fT5|) . i.e., it includes at least the 2 slowest 
non-hydrodynamic modes. 



IV. NON-HYDRODYNAMIC POLES IN AN TV = 4 SYM PLASMA 



The analytical properties of thermal retarded correlators at strong couplin g an d their calculation via the AdS/CFT 
correspondence [l6[ have been discussed in the literature in great detail [2ll |22|. Through the duality [23| . the poles 
of the retarded thermal 2-point function of T xy correspond to the quasinormal frequencies in the so-called scalar 
channel, which amounts to solving the equation of motion for a massless scalar field minimally coupled to gravity in 
the bulk [24[ ■ It was found that for strongly coupled A/" = 4 SYM theory the poles always come in pairs with the same 
imaginary part and o ppo site real parts [25| , which is very different than the weak coupling behavior inferred from the 
Boltzmann equation [IJJ (see Fig. 1). It would be interesting to check if this retarded correlator in weakly coupled 
Af = 4 SYM indeed possesses poles only on the imaginary axis. As shown in (23|, due to rotational invariance, the 
three independent Green's functions that describe the fluctuations become identical at zero wavenumber and, thus, 
in this limit the non-hydrodynamical poles in all of these channels become the same. 

In the supergravity approximation, the retarded correlator of the glueball operator TrF 2 in strongly coupled J\f = 4 
SYM is also found by solving the same equation of motion for a minimally coupled scalar in the bulk. This occurs 
because both operators have their UV scaling dimension equals to 4, which corresponds to a massless scalar field in 
the bulk [l|| . The poles of this glueball correlator determine the masses and the decay properties of the glueballs at 
finite T. It is clear in this context, however, that the fundamental mode of this glueball correlator must indeed be 
doubly degenerate. This occurs because the mass of the state, which is nonzero, only appears as m 2 (i.e, the poles have 
opposite real parts) and in the deconfined phase glueballs must decay (hence the poles must have a nonzero imaginary 
part). Therefore, while the numerical value of the poles of the TrF 2 (or T xy ) correlator will change according to 
the conformal theory used in the calculation, the general structure discussed above concerning the distribution of the 
poles in the complex plane should remain valid for conformal plasmas. Note, however, that while in the supergravity 
approximation the singularities of the Green's functions do not depend on the t'Hooft coupling or N c , the analytical 
properties of the Green's functions are expected to change outside the supergravity limit (see, for instance, the 
discussion in (26|). It would be interesting to investigate the analytical structure of this correlator in non-conformal 
plasmas described by bottom-up gravity theories constructed to mimic some properties displayed by QCD at nonzero 
temperature [27| . 

A crucial insight derived in flOfl, which will be used extensively here, was that the zero wavenumber limit of the 
poles of the Green's function determines the coefficients of the linear terms in the equations of motion of rclativistic 
dissipative fluid dynamics. 



V. LINEARIZED FLUID DYNAMIC EQUATIONS FOR THE M = 4 SYM PLASMA 

We choose to derive the macroscopic equations of motion for W*" via linear response to small metric perturbations, 
as discussed in [2lj |. The linear transport coefficients are determined from perturbations h/ a> of the metric tensor 
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a) 



Imco 



b) 



Imco 



Rcco 



Rcco 



FIG. 1. Analytic structure of the retarded Green's function associated with shear stress in weakly coupled theories (a), based 
on the Boltzmann equation, and strongly coupled theories (b), based on the AdS/CFT duality. 



g^v = rqiiv _|- j n g au ge theory. While this method can be equally used in weak and strongly coupled gauge 
theories, we will focus on the results for strongly-coupled H = 4 SYM. Within linear response, the variation of the 
energy-momentum tensor T^ iV due to the metric perturbations is 



1 f°° 

5T „v (X) = ~ d 4 X' G% a0 (X - X') h aP {X') , 

2 J-oo 



(19) 



where G^ a ^ (X — X') is the retarded Green's function, whose properties will be obtained from the AdS/CFT cor- 
respondence. We consider only the following metric perturbation h xy ~ h xy (t, z), with all other components of the 
metric tensor left unperturbed . In this case, all the other components of 8T^ V decouple from the xy component 
and we obtain the following expression for ST xy 



5T xy (t,z) 



dt' dz' G xyxv {t -t'-z- z') h xy (t', z') . 



(20) 



The energy-momentum tensor T^ u is assumed to have the traditional fluid-dynamical structure shown in Eq. ([T}, 
which then implies that 



ST xy = T xy^u + U ^ _ rpxy^u^ = _p Q h xy + g^a 



(21) 



where 5ir xy is the xy component of the shear stress tensor created by the metric perturbations and Po is the pressure of 
the unperturbed state. We set the shear stress tensor of the unperturbed state to zero. Using the energy-momentum 
equations of motion we arrive at the following equation [l(| 



Sir" = P h xy + / dt' dz' G x R yxv (t -t'-z- z') h xy (t', z') . 

J — oo 

or, equivalently, in Fourier space, Sn xy (t, z) — [l/(27r) 2 ] Jf° duj dke^ l ^ t+lkz Stt xv (uj, k) where 

57T xv (uj, k) = G R {u, k) h xy {u, k) , 



(22) 



with Gr(uj, k) = — Po + G X ft XV {ui, k). Note that G_r(w, k) has the same analytic structure as G x ^ xv {uj, k) because Po 
does not depend on oj and k. 

Given that Gn(u),k) is a meromorphic function with only non-hydrodynamic poles, in the case of homogeneous 
relaxation (where k = 0) one can formally write fiol ] 



-. oo 

G r (lj) = F(u) + — 



(23) 



where F(ui) and / ra (w) are analytical functions (and we used that G* r (lu) = Gr(—uj*)). Performing the Fourier 
transform and picking up the residues one finds (using that h*{uj) = h(—uj*)) 



6***(t) = P h xy (t) + 9(t) J2 |/«K(0)) K v {u n {$))\ e" r "* cos(fi„t + 5 n ), 



(24) 
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where S n is a constant phase shift. Clearly, one must be careful when dealing with the representation above because, 
in general, the sum in the Eq. (|23p may not converge. Since there is only a few explicit examples where all the poles 
and residues are known analytically, in general, the convergence properties of the sum employed in Gr are not known. 
However, note that in the equation for 8it xv (t) derived above, h enters in the coefficients of the sum. Therefore, the 
sum in Eq. (|24p may converge as long as the metric disturbance varies sufficiently slowly in time, i.e., h goes to zero 
fast enough for u ^ 0. Note, however, that this is indeed the case we are interested in since we want to study the 
response of the system to an external agent (the metric variations) that varies sufficiently slow in time (near the fluid 
regime). Thus, it is simple to show that under these conditions, only the first two poles (n — 1 in the sum, i.e, two 
distinct time scales) will contribute at sufficiently long (though finite) times tTi 1. 

Including only the 2 slowest non-hydrodynamic modes close to the origin in the w-complex plane, we obtain the 
following linearized equation of motion for 8ir xv in an AdS/CFT configuration (see Fig. 1-b) directly from Eq. (1221) . 

[$ 2 (0)d t 2 + $!(0)9 t + 1] 5*** = Cx{Q)h xy - ( d " C f k) ) d 2 z h(t, z) 

\ Z J k—0 

+ [C 2 (0) + Ci(0)$i(0)] h(t, z) + 0(h(t, z), dfh(t, z)), (25) 



where we defined 



p 



P 



C p (k) = -dZG R (co,k) , $ x (fc) = — --, $ a (fc) = ... ... . (26) 



=0 Wi(fc) U>2{k) U)i(k)u)2.{k) 



Note that, due to the symmetries of the retarded Green's function, it was not possible to include only one non- 
hydrodynamic mode, as was possible in the Boltzmann equation and happended in Israel-Stewart theory, since the 
first two poles are symmetric relative to the u imaginary axis and are equally distant from the origin in the complex 
w-plane. Eq. (|25[) came directly from the underlying microscopic theory since our starting point was Eq. (1221) . We 
assumed in the derivation of Eq. (|25[) that the Green's function contains only non-hydrodynamic poles (that are 
functions of k 2 due to rotational invariance), Cq(0) = Gr(0, 0) = 0, and we also limited ourselves to only display the 
terms containing at most 2 derivatives. 

It is easy to show that the macroscopic theory in Eq. (1161) when linearized via metric perturbations becomes (note 
that due to these metric perturbations the shear tensor becomes a xy ~ dth xy ) 

[X2d 2 t + xidt + 1] 6ir x y(t, z) = V d t h(t, z) + nc 2 d 2 h(t, z) + V ( Cl + c 2 )d 2 h(t, z) . (27) 

Note that the coefficient £ in Eq. (fTB|) does not appear in this case because of the specific way we chose to disturb the 
metric. However, it can be shown by calculating the general dispersion relations for the sound and shear channels for 
a conformal fluid that a term like D^D\tt uS > x cannot appear and therefore we must take £ = [28|. 

We can now match the long distance, long time limit of the microscopic theory in Eqs. (|25[) to the macroscopic 
theory in (|2T|) to derive the well-known Kubo formula r\ = Ci(0) = i<9 w G;?(0) and, thus, we recover that rj/s = 

l/(47r) Hi] in the supergravity limit. Moreover, xi,2 = $1,2(0), 2?/c 2 = - (<9 2 C (fc)) = - (<9 2 Gfl(0, k)) 



fc=0 
■xyxy 



and 

k=0 



2r)(c\ — xi) = {d 2 — 9 2 )Gfl(w, k) . Using the results for the Taylor expansion of G x ^ xy derived by [ll| and the 

U) — k — r 1 

calculation of the poles at zero wavenumber from [25], one obtains the following values for the transport coefficients 
in TV = 4 SYM, X i ~ 0.63/(2ttT) and X 2 ~ 0.23/(2ttT) 2 , c 2 = 1/(2ttT) and c x = % x - (2 - ln2)/(27rT). Thus, all 
the coefficients associated with the linear terms in the transient theory in Eq. (|16|) have been determined. There are, 
however, still 10 coefficients in the transient theory that remain to be computed: e^'s, /i's, C3 and C4. Since they 
correspond to nonlinear terms, they cannot be determined from linear response theory. 

As was mentioned before, since the transient theory derived in this paper automatically reduces to the asymptotic 
theory derived in (llj | several results derived within that theory are contained within the general transient theory 
displayed in (fTB|) . For instance, results derived within the fluid-gravity correspondence [2!| or the Bjorken expanding 
systems studied in 30] can be readily recovered. In fact, in the case of a Bjorken expanding system, the transient 
theory defined by Eq. (fT6|) should give the same expressions obtained from the Burnett-like theory in Eq. |9]) at 
sufficiently large times. 

VI. FINAL COMMENTS 

In summary, in this paper we derived the most general equation of motion (see Eq. (116j) ) compatible with conformal 
invariance satisfied by the shear stress tensor of strongly coupled M = 4 SYM theory in the transient regime at 
0(Re~ 2 , Kn 2 , Re _1 Kn 2 ). This equation contains 17 transport coefficients (of which 7 were determined in this paper) 
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and it describes the transient regime experienced by the fluid as it evolves towards its universal asymptotic solution 
given by the gradient expansion computed to C(Kn 2 ) in [ll[- Equation (Tip) is a second-order differential equation 
with respect to propertime for 7T MI/ and, as such, it is structurally different than the relaxation-type equations expected 
to describe the transient fluid dynamics of weakly-coupled systems [Tof ■ 

An important point concerns the stability of Eq. (|16l) with respect to hydrostatic equilibrium. As mentioned above, 
the stability of transient, relativistic fluid dynamics is a nontrivial problem and, in fact, so far only the stability 
conditions for relaxation- type equations have been checked [3] . The inclusion of additional time derivatives affects the 
previous studies and, thus, one must generalize these calculations to verify under which conditions Eq. (|16|) describe 
a stable and causal fluid. 

The novel transient physics contained in Eq. (|16l) may bring some light into the description of the early time 
dynamics of the strongly coupled QGP formed in ultrarelativistic heavy ion collisions. The authors thank H. Niemi, 
H. Warringa, G. Torrieri, and D. Rischke for discussions. We thank the Hclmholtz International Center for FAIR within 
the framework of the LOEWE program for support. J. N. thanks Conselho Nacional de Desenvolvimento Cientifico 
e Tecnologico (CNPq) and Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP) for support. 
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